Suzaku observes weak flares from IGR J17391-3021 
representing a common low-activity state in this SFXT 
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ABSTRACT 

We present an analysis of a 37-ks observation of the supergiant fast X-ray transient (SFXT) 
IGR J17391-3021 (=XTE J1739-302) gathered with Suzaku. The source evolved from quiescence to 
a low-activity level culminating in three weak flares lasting ~3 ks each in which the peak luminosity 
is only a factor of 5 times that of the pre-flare luminosity. The minimum observed luminosity was 
1.3 X 10^^ ergs"'(^^/2.7kpc)- in the 0.5-lOkeV range. The weak flares are accompanied by significant 
changes in the spectral parameters including a column density (A^h = (4.1^0 5) x 10~'^cm"^) that is ~2-9 
times the absorption measured during quiescence. Accretion of obscuring clumps of stellar wind mate- 
rial can explain both the small flares and the increase in A^h- Placing this observation in the context of 
the recent Swift monitoring campaign, we find that weak-flaring episodes, or at least epochs of enhanced 
activity just above the quiescent level but well below the moderately bright or high-luminosity outbursts, 
represent more than 60+5% of all observations in the 0.5-lOkeV energy range making this the most 
common state in the emission behavior of IGR J 1739 1-3021. 

Subject headings: accretion, accretion disks ; gamma-ray J: general ; stars: neutron ; supergiants ; X-rays: binaries ; 
X-rays: individual (IGR J17391-3021 = XTE J 1739-302) 



1. Introduction 



Su p ergiant Fast X - ray Tr a nsients ( SFXTsi l Smith et al 
20041: ISguera et al.l l2005l l2006t iNegueruela et alJ 
2006bl) form a new subclass within the group of high- 
mass X-ray binaries ( HMXBs) Prev iously, HMXBs 
were catalogued (e.g. iLiu et al. I l200d and references 
therein) into one of either two groups: the tran- 
sient sources with B-emission line stellar compan- 
ions (BeXBs); or the persistent, but variable, HMXBs 
with supergiant O or B-type companions (SGXBs). 
The SFXTs share properties from both groups: like 
BeXBs, they are characterized by outbursts of short 
(hours) to long (days) duration whereby the peak in- 
tensity is several orders of magnitude greater than dur- 
ing quiescence; and like SGXBs, they are systems 
in which a compact object (usually a neutron star) 
is paired with a supergiant OB star Their study is 
therefore important because SFXTs could represent an 
evolutionary link betw e en BeXBs and SGXBs (e.g. 
Negueruela et al.l bOOSi: IChatvl l20ld iBodaghee etal 
2010l:lLiu et al.ll2010l) . 



XTEJ1739-302 ( ISmith et alj 119981) was the first 
hard X-ray transient to display behavior that was 
identified as characteris tic of the class of SFXTs 
dNegueruela et al.l2006bh . It was listed as AX J 1739. 1 -3020 
in the ASCA catal og of X-ray sources in the Galac- 
tic Center region dSakano et alj 120021) . Five years 
after its discovery, INTEGRAL detected a transient 
source named IGR J 1739 1-3021 which, because of 
the positional coincidence and similar emission be- 
havior, was shown to be the har d X-ray (18-50keV) 



counterpart to the RXTE source dSunvaev et alj|2003 



Smith et alJl2003h . A refined X-ray posi tion was ob- 
tained with Chandra dSmith et alJ l2006l) . which led 
to the identification of USNO-B 1.0 0596-0585865 
(= 2MASSJ17391 155-3020380), a supergiant star 
of spectral type 08Iab lo cated ~2.7kpc away, as 
the optical/IR counterpart dNegueruela et aL l2006at 
Chatv et al.ll2008HRahoui et al.ll2008l) . 



This source has been detected numerous times 



ing campaigns with INTEGRAL ( Lutovinov et al. 


2005 


Sguera et al.ll2005t Tiirier et al.ll2007: Blav et al. 


2008 


ChenevezetalJ I2OO8: Romano etal. 2009bt 


Ducci et a 


iJ l2010h and with Swift ( 


Romano et al. 


2008. 


3bt 


Sidolietal. 2008, 2009a b: 


Romano et al. 


2009i 


3, 2010). Most of the outbursts are short (~0.5h) 



and weak (~10cps in 20-40 keV as seen with IN- 
TEGRAL), while a few flares last longer than an 



hour (sometimes days as seen in the 0.5-lOkeV 
band with Swift) and are rather intense (~60 cp s 



or ~400mCrab with INTEGRAL) (iBlay et al.ll2008h . 



The photoelectric absorption is va riable from out 



burst to outburst dSmith et al 
(13!^) X 10 



2006h . reaching up to 



cm ^ during the rising portion of certain 



flaring episodes ( Sidoli et al. 2009allbl) . The dynamic 
rang e between quiescence and peak outburst is ~ 10"^ 
(lin't Z and 20051). Based on multiple observations of 
IGR J1739 1-3021 gathered over the course of a year 



with Swift, iRomano et al.l (l2009al) showed that bright 



outbursts account for only 3-5% of the entire observa- 
tion set, while the source spends 39+5% of the time in 
an inactive state near quiescence. Quiescence, i.e., a 
flux in the 0.5-lOkeV range below the Swift detection 
limit of 1.6 X 10 ' ^ ergcm~^ s~' for a ~1 ks exposure, 
is a very rare state dRomano et al ]l2009ah . 



What triggers the outbursts is still an open ques- 
tion. Several mechanisms have been proposed to 
explain the peculiar X-ray emission behavior of 
SFXTs. The outbursts could be due to one or a 
combination of the following effects: the intermit- 
tent accretion of clumpy wind material by a com- 
pact object in a larger orbital radius than in clas- 
sical SGXBs; passage of the compact object in an 
eccentric orbit through anisotropic winds emanat- 
ing from some supergiant stars; and the dynamics of 
wind accretion onto highly-magnetized neutron stars 



disks, and Rayleigh-Taylor instabilities (e.g.iin't Zand! 


20051 Walter & Zurita Heras 2007: 


Sidohet al.l2007l 


Negueruela et al.ll2008: Bozzo et al. 


2008: Ducci et al. 



i2009, and references therein). 

While the identity of the stellar companion is 
known, the nature of the compact object has yet to 
be conclusively determined. The spectral properties of 
IGR J17391-3021, i.e. its hard power-la w shape and 



possib le spectral cutoff around 13keV (ISidoli et al. 
2009ah . suggest a neutron star primary. However, co- 



herent pulsations or cyclotron absorption lines that 
would firmly establish the presenc e of a neutron star 



have not been detected thus far ( Smith et al.l 11998 



Sguera et al.ll2005l: iBlav et al.ll2008h . Recently, the or- 



bital period of the system was shown to be 5 1 .47(2) d 
with hard X-ray outbursts (>20keV) and quiescent 
epochs coexisting at periastron su ggesting clumpy, 
anisotropic winds ( Drave et"al]|20 1 Ol) . 

In February, 2008, Suzaku observed IGR J 1739 1 -302 1 
for 37 ks of effective exposure time. The data from this 
observation are described in Section |2] Results from 
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Fig. 1. — Image of IGR J 1739 1-3021 as captured by the XISl detector on Suzaku. Pixels have been grouped in blocks 
of 4 (bin: group4 in ds 9). The circular source and square background selection regions are also shown. Coordinates 
are given as equatorial R.A. and Dec. (J2000) where North is up and East is left. 



timing and spectral analyses are reported in Sections 
|3] and |4] respectively. We discuss our findings and 
present our conclusions in Sections |5]-|6] 

2. Observations 

The Suzaku space telescope features two co-aligned 
instruments working in tandem to observe the hard X- 
ray sky ( 0.2-600 keV): the X -ray Imaging Spectrome- 
ter rXIS: lKovama et al.ll2007h and th e Hard X-ray De- 
tector (HXD: iTakahashi et al.ll2007l) . Currently, XIS 
consists of two front-illuminated CCD cameras (XISO 
and X1S3), and one that is back-illuminated (XlSl). 

Observation ID 402066010 (PI: J.A. Tomsick) tar- 
geted 1GRJ17391-3021 on February 22-23, 2008 
(MJD 54518.495-54519.396), for 72.271 ks. After 
removing bad events during data reduction (see be- 
low), and observation gaps due to the occultation of 
the satellite by the Earth, the net exposure time that 
remained was 36.546 ks. The XlS-nominal pointing 
mode was used with the quarter-window option in or- 
der to gain a higher time resolution. 

Data reduction employed v.6.9 of the HEAsoft soft- 
ware package following the procedures described in 
the Suzaku ABC Guide v.2. The task xispi converted 
the unfiltered XIS event files to pulse-invariant chan- 
nels, which were then input into xisrepro to re- 
turn cleaned event files. We extracted the source spec- 
trum from each XIS detector using a circular region of 
120"-radius centered on the position given by Chan- 
dra (see Fig.[T]i. The background was extracted from 
a pair of 240"-wide squares positioned on a source- 



free region of the detector. Based on these source and 
background extraction regions, we generated response 
matrices with xisrmfgen and xissimarf gen. 

To improve photon statistics, we merged the source 
spectral files for the front-illuminated CCDs (XISO and 
XIS3), and we did the same for the background. Each 
bin in the XIS spectra contained a minimum of 150 
counts. Source and background light curves were gen- 
erated in 0.5-lOkeV (XIS) with a time resolution of 
2 s. Timing data from the 3 XIS detectors were com- 
bined into single light curves for the source and for the 
background. 

A tuned (v.2.0) background file and the common 
GTI were used during spectral extraction of HXD- 
PIN. Dead-time corrections were performed with 
hxddtcor. The cosmic X-ray background was mod- 
eled using the flat-field response file from January 29, 
2008, and its contribution was accounted for in the PIN 
background spectrum. Since IGR J 1739 1-3021 was 
very weak during this observation, PIN recorded few 
high-energy photons. Furthermore, IGR J17391-3021 
is located close to the Galactic Center so diffuse emis- 
sion from the Galactic Ridge can influen ce the high 



energ y spectrum of IGR J17391-3021 (lYuasa et al. 



20081). Applying t he best-fitting parameters for the 



ridge emission from lValinia & Marshalll (119981) to the 
PIN spectrum, we find that the ridge contributes over 
50% of the PIN photon counts. Hence, the statistics 
being insufficient for timing and spectral analyses, ob- 
servations with PIN are not discussed any further 
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Table 1: The different states in the emission behavior of IGRJl 7391 -3021 as used in the text. 



State 


absorbed flux {0.5-lOkeV) 
(erg cm~~ s~' ) 


example time 


reference 








(MJD) 








quiescence 


< (1 -2) X 10"'- 


54518.495-54518.920 


this work, see 


also Bozzo et al. 


(2010) 


low 


(2- 10) X 10"'- 


54518.920-54519.396 


this work, see 


also Sidoli e^L 


(2008); Romano et al. (2009a. 2010); Bozzo et al. (2010) 


medium 


- {2 - 100) X 10"^' 


54635 


Romano et al. 


(2009a, 2010) 


high 


> 10-9 


54692 


Romano et al. 


(2009a. 2010) 





full observation 
0,5-10 keV 



k li 



1- 



hard band 
2-10 keV 



i# A *A 



soft band 
0.5-2 keV 




0.5-10 keV 



ZxlO* axlC 4x10* 5x10* 
lime (s) 



xlO* 7x10' 



Fig. 2.— Light curve of IGR J17391-3021 from XIS 
(0.5-lOkeV). The top panel displays the background- 
subtracted light curve from the source, and the bot- 
tom panel shows the background count rate during 
the observation. Also presented are the source light 
curves from soft (S) 0.5-2 keV and hard (H) 2-10 keV 
bands, as well as the hardness ratio which is defined 
as (H-S)/(Hh-S). Each data point collects 240 s worth 
of exposure. After 33 ks of quiescent emission, a 
low state begins at approximately MJD 545 18.92 (in- 
dicated by the vertical dashed line). 

3. Timing Analysis 

Figure|2]presents the light curve of IGR J 1739 1-3021 
in the 0.5-10 keV band from XIS. In this energy range, 
the average source count rate (combining XIS 0, 1, 
and 3) for the entire observation is 0.251 +0.005 cps. 
The best-fitting model to the (absorbed) spectrum of 
the time-averaged Suzaku observation (see Sect.lH 
yields a coefficient of 2.0 x 10"" erg cm"- ct"' which 
we can use to convert instrumental count rates to 
units of 10 " ergcm"^ s Thus, the average flux 



is (5.0 + 0.1) X lO^'^ergcm-^s-i. 

Based on Fig.|2] we can identify two distinct emis- 
sion states from IGR J 1739 1-3021 in this observation. 
During the initial 33ks, IGR J17391-3021 is in a qui- 
escent state in which the average source count rate 
is 0.100+0.009 cps ((2.0 + 0.2) x lO^'^ergcm^^s-i). 
As with previous observations of IGRJ17391-3021 
when it was at or near quiescence ( Sidoh et aP 2008 ; 



Romano et akl l2009ai l2010t IBozzo et al 



2010l) . our 

Suzaku observation shows that even in this state, 
IGR J17391-3021 varies on time scales of a few hun- 
dred seconds. 

Around 33 ks after the start of the observation (cor- 
responding to MJD 545 18.92), the source enters a low 
state in which the mean count rate of 0.358+0.008 cps 
((7.2 ± 0.2) X lO^'^erg cm"- s"') is a factor 3-4 that of 
the flux in the quiescent state that preceded it. The ob- 
servations prior to and after MJD 54518.92 are hence- 
forth referred to in this analysis as "quiescent" and 
"low" states, respectively (see Table[T]i. In addition, 
three flares lasting ~2 ks each are revealed in the light 
curve at f ~ 40, 54, and 70 ks. The peak flux of 
1.3±0.2cps ((2.6 + 0.4) x 10"" ergcm"^ s"') is an or- 
der of magnitude larger than the flux during the quies- 
cent state, but the maxima are still much lower than the 
typical high-state outbursts expected and seen in such 
objects. We note that the peaks are nearly equidistant 
in time being separated by intervals of approximately 
15 ks. 

The light curve was divided into two pieces corre- 
sponding to a soft (S) band between 0.5 and 2keV, 
and a hard band (H) with energies between 2 and 
lOkeV. The hardness ratio can then be defined as 
(H-S)/(H+S) where a negative value indicates a soft 
spectrum, and where a positive value denotes a hard 
spectrum. Figure[2]displays no obvious correlation be- 
tween the hardness ratio and the source intensity. Bin- 
ning the light curve into longer time segments, and 
plotting the hardness ratio versus the source intensity 
shows only a marginal correlation (but within the er- 
rors) between the two parameters at large intensities. 

As part of a monitoring campaign focusing on 
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Fig. 3.— Light curve of IGR J17391-3021 from Suzaku-XlS, Swift-XKT, and XMM-Newton. Fluxes are in units of 
10"" ergcm"^ s"' as observed (i.e. absorbed) in the 0.5-lOkeV energy band, and time is given in MJD (-54000). 
The top panel displays the background-subtracted light curves from the source during Obs. ID 00030987012 from 
XRT (blue circles, 100-s binning) and from our entire XIS observation (red crosses, 240-s binning). Note the half- 
day gap between the observat ions. In the bottom pa nel, the source light curve based on ~ 1 ks observations from the 
XRT monitoring campaign of lRomano et al. I (l2009ah is shown (blue circles) along with their 3cr upper limits for non- 
detections (green downward arrows). The blue upwards arrows designate the maximum flux of high-state outbursts 
detected by XRT that are situated beyond the scale of the graph (from left to right, their values are 125, 12, and 551). 
The average flux from our 37-ks XIS observation is plotted as a single red cross at the intersection of the dotted lines. 
The purple triangle at MJD 54740 corresponds to the average flux in a 31-ks observation with XMM-Newton from 
BozzoetalJ(l20ia) . 



5 



SFXTs, Swift targeted IGR J 1739 1-3021 at various 
times throughout 2008 (iRomano et al.ll2009al) . In par- 
ticular, there is a ~900-s long XRT pointing (Obs. ID: 
00030987012) several hours prior to our Suzaku obser- 
vation. The 0.5-10 keV Swift-XRT detections and 3cr 
upper limits are shown along with the average count 
rate recorded by XIS in Fig.|3] The instrumental con- 
version coefficient for XRT is 1 .3 x 10 '° erg cm"^ ct"' 



dRomano et al.ll2009ah 



4. Spectral Analysis 

The time-averaged spectrum of IGR J 1739 1-3021 
gathered in the 0.5-lOkeV range by XIS is presented 
in Fig.|4] We fit an absorbed power law, employing 

the abundances of Wilms et a l. '("2000^ and the photo- 

ionization cross sections of Balucinska-Church & McCammon 



(Il992h . to the XIS spectra which yields F =1.4+0.1 
and A^H =(3.6+0.4)xl0^^cm"^ (errors are provided at 
the 90% confidence le vel). These are com patible with 
the values derived by ISidoh et al.l (120081) for out-of- 
outburst epochs in a similar energy range with Swift- 
XRT, with an average luminosity that is slightly higher 
at 4.8 X 10^^ erg s~' compared with 3.9 x 10^^^ erg s ' for 
Sidoli et al. ( 2008|). Both studies assu me a source dis- 
tance of 2.7 kpc dRahoui et al. 2008). The model fits 
the data well (xl = 0.82 for 66 de grees of freedom: 
dof). 

Residuals remain at low energies (<lkeV) which 
could hint at the presence of soft excess emission. 
We know that this source exhibits a soft excess 
since this emission was detected with XMM-Newton 
( Bozzo et al. 2010l) . To account for this feature, we 
added a blackbody whose temperature was fixed to the 
best-fitting value from XMM-A^ewfo« (^r = 1.2 keV). 
This reduces the;ify to 0.55 for 65 dof (F-test probabil- 
ity of 5x10-5), and A^h = x lO 
the photon index is poorly constrained (F 

We split the spectra into two segments correspond- 
ing to the quiescent and low states described in Sect.|3] 
(i.e. before and after MJD 54518.92). Table|2l which 
lists the spectral parameters of IGR J 1739 1-3021 for 
each epoch, shows that there is a clear distinction in 
column density, photon index, and luminosity. The A^h 
is at least twice as high in the low state ((4.1^^^) x 
10^^ cm"^) as it is in the quiescent phase ((l.O;^^^) x 
10^^ cm"^). The source spectrum is softer (F = 1.5 + 
0.1) during the low state compared to the quiescent 
state (F = 1.0 + 0.3). The unabsorbed 0.5-lOkeV lu- 
minosity rises from 1.3 when the source is quiescent to 



22 -2 

cm 

)+0.3\ 



but 



: 0.2!yj). 



7.4X10"'-' erg s ' when it is in the low state. 

Figure|5] illustrates the incompatibility (at >99% 
confidence) in the absorbing column from quiescence 
to the low state. It is clear from this figure that the col- 
umn density is the main parameter that evolves from 
quiescent to low states. Indeed, for a constant F ~ 1.3, 
the A^H ranges from ~1 to ~4 x 10^^ cm"^ between qui- 
escent and low states, res pectively. This is different 
from what was seen bv iBozzo et al. ( 2010l) who found 
that the column density remained steady (within the 
statistical errors) during the low state while the F var- 
ied by around 50%. 

To test whether we overestimated the column den- 
sity during the low state by not considering a soft 
excess, we included a blackbody component in the 
spectral fits to quiescent and low states, and we fixed 
its tem perature to the optimal value from IBozzo et"aD 
(1201 oh (1.2keV). The absorbing columns remain in- 
consistent (at the 90% confidence level) between qui- 
escence (A^H < 1 .35 X 10-^ cm"^) and the low state (A^h 
= (1.99!°™) X 10^2 cm-2). When the blackbody is 
fix ed to the lower tem perature of 0.13keV provided 
by Bozzo et al. ( 2010h . the column densities are com- 
patible within the large error bars between the quies 
cent (A^H = (2.5+5-^ 



(5.8! 



-2:0) X 10 



22 cm ^) and low states (A^h 



^^[■|) X 10^"^ cm"'^), but in this case, the result- 
ing blackbody radius is > 100 km which is much larger 
than the expected size of the compact accretor 

For ease of comparison with previous studies, we fit 
the source spectru m with mode l s used for this source 



by other authors ( Sidoli et al.l l2008t iRomano et al 



2009at iBozzo et al.1 120101) . An absorbed blackbody 
provides a better fit than the absorbed power law 
(xl - 0.71 for 66 dof). The parameters from this 
model agree with those describing the blackbody- 
modeled spectrum of IGRJ17391- 3021 when it was 
in the low state as se en with Swift ( ISidoli et al.ll2008 



Romano et al. 2009 al) . The results of this model fit 



to the various states (i.e. time-averaged observation, 
quiescence, and low state) are listed in Table|2] Once 
again, statistically significant differences emerge for 
the column density (as well as for the radius of the 
blackbody emission region) for spectra collected in 
quiescence (A^h ^ 0.22 x 10^^ cm"^) and during the 
low state (A^H 



(1.4!°;^) X 10-2 



cm ^) 



In the low state, a power law with a partial-covering 
absorber and an additional component to account 
for Galactic absorption J^fixed to 1.3 x 10^^ cm"^ 

(5.1!!:°) X 



from iKalberla et all (120051) ) yields A^h 



cm with a covering fraction of 91 + 4% and 
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Fig. 4. — Background-corrected spectrum of 
IGR J17391-3021 fit with an absorbed power law. The 
data represent photon counts from XISl (red), and 
XISO combined with XIS3 (blue). Each bin collects 
a minimum of 150 counts. 

r = 1.76!°-^^ Cr?/dof = 0.56/53). In other words, dur- 
ing the low state that we observed, a large percentage 
of the X-ray emission is shielded by an absorber whose 
A^H is higher than the expected interstellar value. We 
could not constrain the value of the partially-absorbing 
column in the quiescent state. 

A cutoff power law model gives a good fit to the 
time-averaged spectrum (xl/dof - 0.57/65). The 
power law index is F = -0.2^!] $ with a cutoff at 



= 3;^jkeV, and the column density is A^h = 
These are consi stent with 
(I2OIOI) . The in- 



^CUt 

(2.1 

the parameters listed in IBozzo et al 



+ 0.3) X 10^2 cm-^ 



dex and cutoff energy are compatible between the qui- 
escent (F = -0.1 Ecut > 2keV) and low (F = 
-0.3+'^-^, ^cut = 3+j keV) states. In this case, the 
upper limit to the column density during quiescence 
(< 1.5 X 10" cnT^) lies within the 90% confidence in- 
terval of A^H in the low state ((2.2+° !5) x 10^^ cm"-). 

The spectra were also fit with the Mewe-Kaastra- 
Liedahl (MEKAL) model employing a dual absorber: 
one that affects the system's soft X-ray emission from 
hot diffuse gas in the supergiant star's wind (A^Hsg); and 
another that impacts the high-energy power-law emis- 
sion from the compact accret or (Nhx)- The gas tem- 
perature was set to 0.15keV ( Bozzo et al. 2010l) . In 
the quiescent state, we obtain A^Hsg ^ 2.2 x 10^- cm"- 
and A^Hx < 2.4 x 10^^ cm"- for F fixed to the optimal 
value of 1.0 from the absorbed power law fit to the 
quiescent spectrum (xl/dof = 0.54/52). For the spec- 
trum from the low state, A^Hsg = (5.8;^2o) x 10"^ cm"^ 



Table 2: Spectral parameters from a power-law (PL), 
a cutoff power law (CP), and a blackbody (BB) model 
fit to XIS 0-3 spectra of IGR J17391-3021. Fluxes 
(10"'^ergcm"2s"') and luminosities (lO-'"^ ergs"') are 
given as unabsorbed in the 0.5-10 keV range for a 
distance of 2.7 kpc. Quiescence and the low state re- 
fer to pre and post-MJD 545 1 8 .92 spectra, respectively 
(see Fig.|2li. The A^h is given in units of 10^^ cm"^, 
the cutoff energies and blackbody temperatures are in 
keV, and the blackbody radii are in meters. Errors are 
quoted at 90% confidence. 



model 


parameter 


time-averaged 


quiescence 


low state 


PL 


r 

F 
L 

xl/dof 


3.6 ±0.4 
1.4±0.1 

5.5 

4.8 
0.82/66 


1 0+°'' 
l.u „5 

1.0 ±0.3 
1.5 
1.3 
0.68/8 


4 l+°-5 
^■'-0.4 

1.5 ±0.1 

8.5 

7.4 
0.73/54 


CP 


Nil 

r 

Ecat 
F 
L 

xllAof 


2.1 ±0.3 

"0-2-a8 

4.0 
3.5 
0.57/65 


< 1.5 
-0.1 ± 1.3 
> 2 
1.3 
1.1 
0.64/7 


f -)+0.7 

"0-3!?:^ 

3!| 
5.7 
4.9 
0.44/53 


BB 


Nil 

«bb 

F 
L 

xllAoi 


1.2 ±0.2 

1 40+0.06 

72 ±5 

3.4 

3.0 
0.71/66 


< 0.22 

1 <:7+0.16 
'•-"-0.15 

37!« 
1.2 
1.0 
0.73/8 


1 40+0.07 

88!] 
5.1 
4.4 
0.51/54 



Table 3: Same as Table|2] but for key parameters from 
rate-selected spectra from epochs of low count rate 
(LCR), medium count rate (MCR), and high count rate 
(HCR). 



model 


parameter 


LCR 


MCR 


HCR 




count late range (cps) 
average count rate (cps) 
exposure time (ks) 


<0.1 
0.05 
22 


0.1-0.3 
0.19 
11 


>0.3 
0.45 
4 


PL 


/Vh 

r 


< 0.9 
0.6!g;^ 
1.2/11 


4 r' ° 

1.7 ±0.3 
0.83/34 


1.4 ±0.2 
0.50/38 


CP 


Nn 
T 
Fc\a 
xl/dof 


< 0.2 
-0.8 ±0.2 
3 (fixed) 
1.2/10 


2 1+'- 

-op 
-0 1+"^" 

3 (fixed) 

0.71/34 


9 c+0.8 

•^■^-0.6 

-0.3 ±0.2 
3 (fixed) 
0.27/38 


BB 


Nh 
kTub 
xl/dof 


< 0.1 
1.6/10 


'■^-04 

I.4±0.1 
0.72/34 


1 9+0.5 
'—-04 

1.4±0.1 
0.72/34 



and A^Hx < 1.4 x lO^^ cm"^ for F = 1.8 + 0.2 (xl/dof 
- 0.55/65). This suggests that not only is there mate- 
rial around the binary system absorbing X-rays, there 
is additional matter located near the compact object. 
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2.0 



1.0 



0.5 



0.0 





low state 



quiescence 



2 4 6 

N„ (10" cm"^) 

Fig. 5. — Parameter space of A^h and F derived 
from absorbed power laws fit to the spectra of 
IGR J17391-3021 during quiescent and low states. 
The best-fitting parameters are represented by crosses 
and the contours denote confidence levels of 68%, 
90%, and 99%. 



We sought to determine how the spectral parame- 
ters varied with the source intensity. So we created 
spectra from time intervals in which the count rates 
(CR) were <0.1 cps ("low CR" or LCR), between 0.1 
and 0.3 cps ("medium CR" or MCR), and >0.3cps 
("high CR" or HCR). Then, each rate-resolved spec- 
trum was fit with an absorbed power law, an absorbed 
cutoff power law with E^ut fixed to our best-fitting 
value (3 keV), and an absorbed blackbody model. Ta- 
ble[3]lists the results of the spectral fits. Regardless of 
the spectral model, the A^h is significantly smaller dur- 
ing the LCR epochs. The MCR and HCR epochs have 
spectra that are consistent, and both are statistically in- 
compatible with the spectrum from LCR epochs. 

5. Discussion 

5.1. Weak Flares and Low Activity 

Our analysis of a Suzaku observation of IGR Jl 7391- 
revealed the source evolving from quiescence to a 
low-activity state. Yet even at this low-flux level, 
the emission from IGR J17391-3021 is characterized 
by variability on time scales of hundreds to thou- 
sands of seconds indicating an accretion that is con- 
tinuous and dynamic. The average luminosity of 
IGR J17391-3021 during the entire Suzaku obser- 
vation is 4.8 X lO^'-'ergs"' (unabsorbed) and 3.1 x 
ergs ' (absorbed) in the 0.5-lOke V range, and 



10^3 

assuming a source distance of 2.7 kpc jRahoui et al 



20081). 

During the quiescent state, which consists of the 
first ~33 ks of the observation, the luminosity in the 
0.5-lOkeV band is only 1.1 x lO-'^ergs"' as ob- 
served, and 1.3 X lO^^ergs"' when corrected for ab- 
sorption. This is e quivalent to the low est luminosity 



detected by ASCA dSakano et al.l 120021) . In our flux 



limited spectrum (LCR), the unabsorbed luminosity 
is 8.7 X lO-'^ergs"', i.e., twic e the minimum va lues 
meas ured wifli XMM-New ton (IBozzo et al.ll2010l) and 
Swift dRomano et al.ll2010l) . 

Concluding our observation is a ~39-ks epoch of 
enhanced activity including three flares lasting 1-3 ks 
each where the peak luminosities are only a factor 
of ~5 greater than the quiescent baseline emission 
that preceded it. Clearly, these flares are significantly 
weaker than the typical high-luminosity outbursts of 
SFXTs in which the dynamic range is a few orders 
of magnitude. These small spikes are reminiscent of 
the short, low-inte nsity flares seen in this source by 
Bozzo et al] ( l2010b with XMM-Newton. 

Separated by regular intervals of approximately 
15 ks, the flares appear to be quasi-periodic. This 
timescale is well short of the estab lished orbital pe- 
riod of 51 days ( Drave et al. 2010l) . and it is prob- 



ably too long for a pulsation per iod if the compac t 



Ducci et all (l2010h 



object is, indeed, a neutron star, 
detected quasi-periodic outbursts of similar frequency 
from this source in INTEGRAL observations, but for 
longer flares (between 5 and 10 ks) of much higher 
intensity (Lx ~ lO-'^ erg s"'). They ascribe the quasi- 
periodic flaring to the formation and dissipation of 
an unstable accretion disk which is the mechanism 



favored by ISmith et al.l (l2006h to explain bright out- 
bursts in this source. The pulsation period of the sus- 
pected neutron star has never been measured, despite 
numerous observations in and out of outburst. If the 
pulse period of the suspected NS in this system were 
known, then a change in the sign of the period deriva- 
3021 tive c ould indicate a di ssipation and reformation of a 
disk jTaam et al.lll988h . This 15-ks frequency could 
possibly correspond to a periodicity but we would 
require further observations to confirm or reject this 
possibility. 

5.2. Spectral Evolution 

The evolution of the emission from quiescence to 
the low state is accompanied by changes in the X-ray 
spectral parameters. The A^h is higher and the spec- 
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tral slope is steeper during the low state than during 
quiescence. Similarly, the high count-rate resolved 
spectrum showed larger A^h and photon index than the 
low count-rate resolved spectrum. The main parame- 
ter that evolves is the column density. This is diff erent 
from the weak flares seen bv lBozzo et al. I (l2010l) who 
found that as the average count rate of the spectrum in- 
creased, r became smaller (i.e. flatter, harder spectra, 
and more positive hardness ratios) while the A^h did 
not change significantly. 

Differences in the spectral properties of the flares 
caught by Suzaku and the ones seen with XMM- 
Newton could be due to unequal geometric configu- 
rations of the accretion stream between the observa- 
tions. The capture by the NS of a clump of wind mate- 
rial will lead to an increase in the mass-accretion rate 
which both instruments detect as a spike in the X-ray 
luminosity. A partial-covering model fit to the X-ray 



spectr a from the weak flares observed by iBozzo et al. 



( 2010l) suggests that up to 40% of the X-ray emission 
is not absorbed by the clump itself, but rather by the 
interstellar medium. The same model fit to our low- 
state spectrum reveals an absorber that shields a much 
larger fraction of the X-ray source (less than 13% of 
the emission is unblocked). 

On the other hand, if prior to being accreted, the 
captured material crosses the line of sight between 
the observer and the NS, the flare will be accom- 
panied by an increase in the column d ensity, which 



is pre cisely what we observe (see also iRampv et al. 



20091) . The expected interstellar value along the line 



of sight is (1.2-1. 3)><10^^ cm~^ dOickev & Lockman 



199d iKalberla et aljl2005h . consistent with what 



measured during quiescence. This represents an up- 
per limit in the direction of the source since it results 
from the integration of absorbing material through the 
entire Galaxy, while IGR J 1739 1-3021 is located only 
~3 kpc away. The A^e from a power law fit to the spec- 
trum of IGRJ17391-3021 during its low state is 2-A 
times the interstellar value and 2-9 times the quies- 
cent value. This suggests that the source is not strongly 
absorbed (intrinsically) during this observation except 
when the NS accretes a clump passing along our line of 
sight. Material in a clump located near the NS will be 
photo-ionized so the true column density is expected 
to be higher than our measured value which assumes 
cold neutral gas. Our largest column density remains 
an order of magnitude below the maxim um A^h mea 
sured for this source, (32 + 3) x 10^^ cm (iSmith et al. 



KXTE that did not extend below 2.5 keV. 

If the origin of the hard X-ray emission is a black- 
body, then this radius expands from 'il^^^m. during 
^^^m when the source is in the low 



quiescence to 

state, in agreement with iRomano et al.l ( l2009al) who 
showed the blackbody radius growing with intensity. It 
is tempting to associate this small blackbody-emitting 
region with the polar caps of a magnetized accret - 



ing neutr on star (^bb ^ O.I/ ^ns^ iHickox et al.l 120041) . 
However, Bozzo et al. ( 2010l) caution that models de- 
scribing spherically-symmetric accretion onto magne- 
tized neutron stars predict an inverse relation between 
the size of the hot spot and the square root of the X- 
ray luminosity, at least for sys tems with an average 
Lx < 10^'' erg s ' dWhite et al.lll983i, and references 
therein). 

5.3. The Ubiquity of Low Activity 

Our Suzaku observation of IGR J17391-3021, 
when placed in the context of the long-te rm monitoring 



campa ign of this source by Swift-XKY jRomano et al, 



2009ai) . shows that the average flux that we recorded 
is consistent with many of the XRT detections at a low 
flux level (F ~ (2- 10) x 10^'^ergcm"^ s"'), and is 
on par with a few of the 3cr upper limits to the count 
rates of non-detections as measured by XRT (Fig.O. 

Weak flares that are just above quiescence, and 
whose peak fluxes are at the level of the faint de- 
tections in the Swift monitoring data, compels us 
to propose four states in the emission behavior of 
IGRJ17391-3021: the extremely low but variable 
emission of quiescence (average Fabs ^ (1 - 2) x 
10"'~ergcm"^ s"' in 0.5-lOkeV); a low-activity state 
that potentially includes we a k flar es such as those 
seen here and in Bozzo et al. I (1201 oh {F ~ {2- 10) X 



10" erg cm"" s" ); a medium state with flares of 
larger amplitude {F ~ {2 - 10) x 10"" erg cm"- s"'); 
and finally, a high state in which the outbursts repre- 
sent the maximum source luminosity {F > 10"'' erg cm"^ s"' ) 
The boundaries are not strict since the peak fluxes 
in the low state extend up to 3 x 10"" erg cm"- s"', 
but they provide a useful classification for the various 
emission states in this source. 

If we assume that the Swift detections when the 
count rates were sufficiently low (i.e. F ~ [2 - 10) x 
10"'- ergcm"^ s"') were due to low states that may or 
may not include weak flares, then we can define a low- 
state duty cycle (DC) as follows; 



2006h . although their result relied on a spectrum from 
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DC- — =60: 



:5% 



(1) 



Here, Srio^ is the sum of exposure times of all ob- 
servations in which IGR J 1739 1-3021 is detected at 
an average flux within the range of F ~ (2 - 10) x 
10~'^ ergcm"^ s"', including the Suzaku and XMM- 
Newton observations, while YTioi denotes the total ex- 
posure from all observations. As expected, the duty 
cycle is not far from the val ue of 61+5% provided in 
Romano et all ( l2009all201ol) . 

Low-activity states in which the flux is above the 
"typical" quiescent level, but well below the medium 
and bright flaring episodes, are therefore quite com- 
mon in this source. In the field of IGR J17391-3021, 
the Umiting flux for XRT is 1.6 x 10 and 1.5 x 
lO''^' ergcm '^ s~' for exposure tim es of 1 and 30 ks, 
respectively jRomano et al.ll2009ah . Therefore, most 
Swift upper limits would be detections if the exposures 
had been equivalent to those of the Suzaku and XMM- 
Newton observations suggesting that the duty cycle of 
the low state is much higher 



According to the source ephemeris jDrave et al. 



20101) . the midpoint of our observation coincides with 
an orbital phase of 0.37+q "j whereas the XMM-Newton 
observation took place at phase 0.8. This is another in- 
dication that the low state is common since it is not 
confined to a specific part of the orbit. 

Despite the large duty cycle, weak flares lasting 
a few hundred to a few thousand seconds can be 
missed by short monitoring observations. Because 
their peak flux is so low, detecting a low state re- 
quires prolonged observations where the long inte- 
gration time guarantees that we will achieve the sen- 
sitivity necessary to detect such low levels of activity. 
Given that IGR J17391-3021 is a prototypical mem- 
ber of its class, it would not be surprising to find weak 
flares from other SFXTs as telescopes continue to de- 
vote observing time to such objects. In fact, recent 
Suzaku and XMM-Newton observations of the SFXT 
IGR J08408 -4503 included simil ar low states that 
were punctuated by weak flares ( Bozzo et al. 2010t 
Sidoli et alJl2010l) . 



6. Summary & Conclusions 

The main results from our Suzaku observation of 
tiie SFXT IGR J 1739 1-3021 can be summarized as 
follows; 



- For most of the observation, IGR J 1739 1-3021 
is in what we believe represents the quiescent 
state. The unabsorbed luminosity in the 0.5- 
lOkeV range is 1.3 x lO"'-' erg s"', and as low 
as 8.7 X 10^*^ ergs"' in the flux-limited spec- 
trum. The spectrum from this epoch is well de- 
scribed by an absorbed, hard power law with A^h 
= (1 .0+JJ-^) x 10^2 cm-" and r = 1 .0 + 0.3, or by 
an absorbed thermal blackbody with a tempera- 
ture of 1.57+|J[^keV. 

- Around 33 ks into the observation (MJD 545 1 8.92) 
the source became slightly more active ("low 
state") culminating in a series of weak flares 
each lasting ~3ks whose peak luminosity is 
only a factor of ~5 greater than the quiescent 
emission that preceded it. During this low state, 
the A^H reached (4.I+Q4) x 10^- cm"- which is a 
factor of ~2-A times that of the expected inter- 
stellar value, and which is ~2-9 times the A^h 
measured during quiescence. The accretion of 
small clumps of stellar wind material can be re- 
sponsible for the weak flares, while the increase 
in A^H can be explained by the passage of the 
clump in the line-of-sight to the X-ray source 
prior to accretion. 

- The average luminosity of the entire observation 
is 4.8 X lO^-' erg s ' . When placed in the full Ught 
cu rve from the ^w i Yit-XRT monitoring campaign 
by iRomano et alj (l2009ah . we find that this lu- 
minosity is on par with the low-intensity XRT 
detections (as well as a few of the 3cr upper lim- 
its from non-detections). Furthermore, we sus- 
pect that many of the faintest detections from the 
XRT contain similar low states (potentially in- 
cluding weak flares) that are just above the qui- 
escent flux. We estimate that the duty cycle of 
the low state is 60+5%, which makes it the most 
common state for this source. 

As a prototypical member of the class of SFXTs, 
IGRJ17391-3021 holds valuable clues to the accre- 
tion processes of these intriguing transients. Contin- 
ued monitoring of this source by sensitive instruments 
such as Suzaku should fuel the ongoing debate con- 
cerning the mechanisms that are responsible for their 
peculiar emission behavior. 
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